A truth about Brownian motion in gases and in general 
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Real thermal motion of gas molecules, free electrons, etc., at long time intervals (much greater 
than mean free-flight time) possesses, contrary to its popular mathematical models, essentially non- 
Gaussian statistics. A simple proof of this statement is suggested basing on only the determinism 
and reversibility of microscopic dynamics and besides incidentally derived virial expansion of a path 
probability distribution of molecular Brownian particle. 
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1. Thermal chaotic Brownian motion of particles of the 
matter is mechanism of diffusion as well as many other 
transport processes. It is interesting by itself too, since 
statistics of random walk of test or marked "Brownian" 
particles contains complete information about transport 
process as a whole including its noise and fluctuations. 
In spite of this, dynamical theory of molecular Brown- 
ian motion remains almost undeveloped. Apparently, the 
common opinion is that anyway it should confirm the well 
known beautiful probability-theoretic scheme based on 
the celebrated "law of large numbers" [l[ which foretells 
that at sufficiently large spatio-temporal scales probabil- 
ity density of displacement, or path, Aii, of a Brownian 
particle (BP) during time t has universal form of the 
Gaussian distribution. For symmetric 3-D random walk 
that is 

V G (t,AR) = {AirDt)' 3 / 2 exp (-AR 2 /4Dt) . 

In fact, it was rigorously proved 0, 0| that such is 
the asymptotic of chaotic walk of hard ball in so non- 
random environment as static periodic lattice of scat- 
terers (also hard balls). Seemingly, Brownian motion 
of gas atoms or molecules all the more can not display 
something else. Indeed, the Boltzmann-Lorentz equation 
[3| , which comes from the Boltzmann equation when ap- 
plied to self-diffusion of gas atoms, produces the Gaussian 
asymptotic. Nevertheless, not all is so simple as seems. 

Recall, firstly, that the Boltzmann-Lorentz equation 
along with the Boltzmann's one is rather rough model 
since results from violent truncation of the exact BBGKY 
hierarchy of equations and its closure with the help of ar- 
chaic Boltzmann's "molecular chaos" hypothesis [j, [H, [|| . 
However, it was pointed out more than once that 
Boltzmann's argumentation 0] looks enough convincing 
only if applied to spatially uniform gas. As to attempts to 
deduce "molecular chaos" from BBGKY equations them- 
selves, they proved to be unsuccessful In opposite, it 
was shown [1(| (see also [ll|, [l2j]) that reformulation of 
BBGKY equations in terms of inter-particle "collisions" 
(instead of continuous interactions) leads to equations 
which are evidently incompatible with this hypothesis if 
a gas is spatially nonuniform. In other words, the very 
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fact of participation of particles in the same collision (or 
connected conglomerate of collisions) at indicated place 
is sufficient reason for mutual statistical dependence of 
the particles. At the same time, if wishing to consider 
displacement of BP (test particle) one has to localize its 
initial position and thereby disturb uniformity (transla- 
tional invariance) of infinitely many distribution func- 
tions (DF) what describe BP under its interactions with 
gas. Consequently, after all one has to deal with an infi- 
nite chain of equations. Attempts to construct approxi- 
mate solutions to them were made in works [lfjj (see also 
(Tl| ) and [l2T |. Their findings, of course, qualitatively 
differ from the Gaussian asymptotic (but qualitatively 
confirm the early phenomenological theory [13l Il4 Il5||). 

Secondly, Gaussian asymptotic of the ball's walk in 
static lattice is closely connected with its ergodicity 
0) S EB] : results of time averaging over n — ► oo frag- 
ments of its trajectory (with each fragment consisting of 
free flight and collision) almost surely are independent on 
the trajectory. This is possible due to the fact that any 
trajectory is fully determined by its single fragment (e.g. 
initial one), and a ratio of number, d — 5n , of quantities 
describing the trajectory to number, p — 5 , of its specify- 
ing parameters is unrestrictedly large, d/p — ► oo . But a 
(weakly non-ideal) gas is quite another matter. Here any 
trajectory of Brownian (test) particle is made of n rj t/r 
similar fragments, with t being observation time and r 
mean free flight time, and is exhaustively characterized 
by d = 6n quantities. However, it is specified by not only 
initial state of the particle under consideration but also 
initial states of many other particles. It is easy to see 12 1 
that a number m of such particles grows with time by far 
faster than n , so that n/m — > . Thus now the ratio of 
number of quantities what completely describe all details 
of a trajectory to number, p = 6m , of parameters (initial 
conditions) what determine these details, tends to zero, 
d/p — ► , as contrasted with above case. By this reason, 
hardly results of time averaging (for instance, "diffusiv- 
ity" or "collision probabilities") will turn out to have 
"almost sure" same limits for different trajectories. Here 
are no grounds for ergodicity, the "law of large numbers" 
(which needs in beforehand prescribed "collision proba- 
bilities" []], [TT|, [13]) an d Gaussian asymptotic. Even in 
finite gas of N < oo particles (in box or torus) necessary 
grounds for ergodicity arise not earlier than after time 
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> Nt [Hj]- This gives illustration of footnote remark 
in [l(| that a value of ergodicity in physics looks over- 
estimated since limit N — ► oo is much more important 
for physics than limit t — > oo . Besides, in essence, such 
warning was highlighted already by N. Kryloy |18l|. 

In principle, similar reasonings jTH, [HI, 0, l2d||. as well 
as BBGKY equations or their analogues, are applicable 
to any realistic many-particle system, and therefore par- 
ticular case of gas has general physical significance. But 
BBGKY equations yet are so bad investigated that none 
supplementary tools would be superfluous. Below, bas- 
ing on determinism and reversibility of microscopic dy- 
namics only, we will derive a kind of virial expansion for 
true probability distributuion of the BP's displacement, 
Vb (t, AR) . Then connect coefficients of this expansion 
with join correlation functions of BP and gas and with 
usual join DF. Finally, show that all these taken together 
imply a differential inequality for Vo(t,AR) which def- 
initely forbids Gaussian asymptotic of V {t, AR) (but 

allows for non-Gaussian distribution obtained in [12| by 
approximative solving BGKY equations). 

2. Let a gas of N 3> 1 atoms in volume il contains a 
"Brownian particle" (BP). Consider, under the thermo- 
dynamical limit N — > oo, Q — > oo, N/Q = v$ = const , 
statistical ensemble of phase trajectories of this system 
which responds to canonical equilibrium distribution of 
its initial state at t — 0, supposing that at t > BP is 
subjected to constant external force / . Time reversibil- 
ity of phase trajectories implies many relations between 
statistical characteristics of their ensemble [2l|, l22|, |23| . 
In particular, according to [211 ] one can write 

(A(q(t))B(q(0)) e £ ^/ T ) = (B(q(t))A(q(0))) (1) 

Here brackets (...) designate ensemble averaging, q(Q) = 
q = {R, ri, ... ,r/v} are space coordinates of BP and 
atoms, S(t) = f ■ [R(t) — R(0)] is work made by the 
external force during time interval t , A(q) and B(q) are 
"arbitrary functions" , and T is initial temperature of the 
system. This equality holds also when BP has internal 
degrees of freedom. 

Let B(q) = n5{R - Ro) exp [- £V U(rj)/T] and 
A(q) = S(R - R') . Then r.h.s. in (TJ takes form 

(B(q(t))A(q)) =F{t,Ro,cl>\R'} = V (t, Ro - R') + 
+ E£=i("i5/»l) I" Mt, Ro, n ...r n \R>) l£ =1 <f>(rj) , 

where <j){r) = exp [— U(r)/T] — 1 , symbol J n means in- 
tegration over n...r n , function Vo{t, Rq — R') is prob- 
ability density of finding BP at t > at point Ro, and 
F n (t, Ro,ri ... r n \R') probability density of this event 
and simultaneously finding some atoms at points rj , 
under condition that BP had started from point R' . 
Clearly, V o (0, r) = S(r) , J V (t, r) dr = 1 , and 

F»(0, Ro, n... r n \R') = 6{Ro - R') F r [ eq) ( ri ... r n \R ) , 
F{0, R Q , <t>\R'} = 5{Ro - R') T^{^\Ro} , 
^{<j>\Ro} = 1 + 
+ Zn=i("oM) S n F^\ ri ...r n \Ro) UU 0(r 3 -) , 



where F„ (n... r n \Ro) are equilibrium DF and 
{4>\Ro} their generating functional. In respect to 
atoms all the DF are normalized as usually [H, @ , that is 
F n (...r k ... \R') -* F n -i(...r k -i,r k+ i... \R') if r k oo , 
and F[ eq \ri\R') — > 1 if v\ — > oo . At one this require- 
ment expresses the "principle of decay of correlations" . 

Initial DF -F„(0, Ro, T\ ... r n \R') include all equilib- 
rium correlations between BP and atoms. But any evo- 
lution of (ensemble of states of) the system is conju- 
gated with disturbance of detailed balance of collisions 
and therefore births additional correlations due to joint 
participation of particles in excess collisions (or "joint 
nonparticipation" in missing ones). These correlations 
could be qualified as non-equilibrium ones, but with those 
reservation that sometimes they describe evolution of not 
so much a system on its own as our information about it. 
Traditionally, their contributions to DF are termed "cor- 
relation functions" (CF) 0, 0, [24| • Let us designate 
them as V n (t,Ro,r\ ...r n |i2') and pick out them from 
F n (t > 0, ...) by definition as follows: 

T{t, Ro, 4>\R'} = T^{<t>\Ro} [ Vo{t, R - Rf) + 
+ ELW»!)f^(^,ri...r I1 | J R')n; =1 ^)] 

From here it is obvious that V^(0,...) = 0, and 
V n (t,...) — > if at least one of the points r k — > oo . 
In particular, 

Fi(*,B ,ri|JI / ) = 
= Vo(t, Ro - R')F[ eq) (n \Ro) + V 1 (t, Ro, ri \R') , 

where function V\ (t, Rq, ri\R') comes from the usual 
"pair correlation function" @, [24j in the full two-particle 
phase space: 

V 1 (t,.Ro ) r 1 |Jl') = / /Vi(i,i?o,ri,P ,Pi|-R')^Podpi 

(2) 

The latter can be connected with the full DF exactly like 
V\{t, Ro, ri\R') is connected with coordinate ones: 

Fi(t,Jlo,ri ) J%,p 1 |Jl') = V {t,Ro-R',Po)x 
x F[ eq) (n , pi \Ro) + Vi (t, Ro, n, P , pi \R!) 

Next, consider left side of ([T]). Multiply and di- 
vide it by (B(q)) and use the fact that any expression 
(<f> B(q)) / ' (B(q)) , with $ being some functional of sys- 
tem's phase trajectory, in view of determinism of transla- 
tion along any concrete trajectory can be interpreted as 
<&'s average over new ensemble of trajectories induced by 
new, non-equilibrium, ensemble of initial coordinates and 
momentums of the system, namely, defined by probabil- 
ity distribution p{q,p) — B[q) po(q,p)/(B(q)) , where 
Po(q,p) denotes original equilibrium distribution. Notic- 
ing also that (B(q)) = .F^-J^l-Ro} , we come to 

(A{ q {t))B{ q {Q))e- £ ^' T ) = 
= V{t,R'\<t>,Ro} e -f-lR'-K>]/T T^{(p\Ro} , 

where V{t, R'\4>, Ro] is probability density of finding 
BP at time t at point R' ander conditions that initially 
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it was placed at point Rq while the gas was in such per- 
turbed spatially nonuniform non-equilibrium state which 
would be equilibrium in presence of the potential U (r) . 
Mean concentration of atoms in such state equals to 

v{r\<f>, R } = [1 + 4>{r)]6\nF (eq] {4>\Ro}l54>{r) 
Thus, in total, we arrive at exact relation 

V{t,R'\cf>,K Q }e-^ R '-^/ T = V Q (t,Ro-R') + 
+ En=i(KM)I n V n (t 7 R ,r 1 ...r n \R')U; =1 Hr J ) 

(4) 

which connects probability distribution of BP's path in 
initially non-equilibrium nonuniform gas and analogous 
distribution, together with generating functional of "non- 
equilibrium correlations" between previous path of BP 
and its current environment, for initially equilibrium uni- 
form gas. In case / = , therefore, r.h.s. represents 
wholly equilibrium Brownian motion. Notice that simi- 
lar relation, for perturbed gas density in place of the BP's 
distribution, was considered in [l7j |. 

3. To understand benefits of this relation, it is useful 
to consider it under the formal "Boltzmann-Grad limit" 
Q when gas density increases, vq — ■> oo , while radii 
of (short-range repulsive) interactions between BP and 
atoms, r& , and inter-atomic, r a , decrease in such way 
that "gas parameters" r^u and r%vo go to zero but 
mean free paths of atoms, A = (7rr^o) _1 , and mean free 
path of BP, A = (tit^o) -1 j stay fixed. Thus one obtains 
"ideal weakly-non-ideal gas" where Fn Q ^ — ► 1 , in those 
sense that, for example, vq /[-F^ \tx\Rq) — l]dr\ ~ 
~ v o r l ~ * ■ As the cosequence, 

;y{r|0,i? o } -> v [l + tftr)] = u exp [-U{r)/T] 

This simplification makes it quite obvious that in the 
Boltzmann-Grad limit on left side of (0J effects of the 
density perturbation depend on its relative measure cp(r) 
only but not on vq . Hence, r.h.s. of (01 also depends on 
4>(r) only, that is at any given function </>(r) (and fixed 
A and A ) there exist finite limits 

lim vl f n V n (t,R Q ,r 1 ...r n \R')U% 1 (f>(r j ) + 0, oo 

Moreover, choosing 4>(r) = 4> = const inside suffi- 
ciently large sphere (for instance, r — Rq\ < v s to with 
v s being sound speed in our gas and to > t ), such that a 
fortiori none correlations go beyond it, we can conclude 
that (at fixed A , A and t ) there exist limits 

lim ^J n V n {t,R ,r 1 ...r n \R) = V n (t, Rq - R') 

(5) 

At that, the gas described by left side of ^ effectively, 
from the point of view of BP, is spatially uniform at start 
and what is more during any given time of keeping BP un- 
der observation. But its density is (1 + 0) times greater 
than density of the right-hand gas. Correspondingly, ^ 
takes the form of simple virial expansion: 



V (t,-AR;l + 4,)ef- AR / T = ( . 
= E,T=o 4> n V n {t,AR;l)/nl W 

Here we introduced AR = Rq — R and third argu- 
ment representing (dimensionless) gas density measured 
in units of v , so that V n (... ; 1) = V n {... ) . 

As it follows from (JH])-©, in weakly nonideal gas spa- 
tial non-equilibrium correlations do not disappear in rela- 
tive sense, i.e. if counting upon (elementary volume v^ 1 
what fall at) one atom, although they disappear in literal 
sense: V n (t, Rq, n... r n \R') — > . 

Let us interpret this statement on the base of those few 
things what are known about the correlation functions 
SSMGlEiHl, first of all, Vi{t,B^Ti,P^,jn\B!). 
In respect to n — Rq this pair CF is nonzero inside 
a "collision cylinder" having radius and directed in 
parallel to relative velocity of the pair of particles (BP 
and some atom) where absolute vlaue of the CF is com- 
parable with product of one-particle DF, i.e. with first 
term on r.h.s. of ([3]). However, the approximative theory 
1Q, 24], conventionally accepted as satisfactory one, says 
nothing about length of the cylinder. On the other hand, 
formula ([5]) clearly prompts that in rigorous theory effec- 
tive length of collision cylinder is finite and has size of 
order of A . Indeed, let it be ciA , then volume of the 
cylinder equals to nrf c\A = c\Uq , and 

fVifaRo^^PoiPilR^dn = 

= c x UQ l Vx{t,AR,P , Pl ) , {<) 

where V\(t, AR, _Fb?Pi) represents mean value of 
Vx{t, Rq, r%, Po,pi|.R') within the cylinder. Combining 
this expression with ^ and j5f, at n = 1, we make 
certain of finitness of the limit: 

Vi(t,AR) = a J JV 1 {t,AR,P ,p 1 )dP dp 1 (8) 

At the same time, in literal sense formulas ^ and ([5]) 
yield Vi{t,Ro,n\Bf) ~ Vt(t, AR) r b 2 /|n - P | 2 ~ 
Vq 1 Vi(t, AR)/(A\ri - R \ 2 ) -> , since under inte- 
gration over momentums a very small part of differently 
oriented collision cylinders only covers (fixed) point r± . 

4. We have come to curious conclusions. 

Operating with coordinate r\ , let us average identity 
(O over the collision cylinder: 

F x {t,AR,P , Pl )_ = 
= V (t,AR,P )G^( Pl ) + Vi(t,AR,P , Pl ) 

Here G lyeq \p\) is the Maxwellian distribution, and the 
result, Fi(t, AR, Pq; Pi) i can be treated as ensemble av- 
erage of density of number of two-particle configurations 
matching up collisions [nil ] . Then, performing integration 
over momentums and applying formula we obtain 

Wi(t,AR) = J /Fi(t,Afl,fl,,pi)dfl,dpi = , . 
= V (t,AR) +c^ 1 V l (t,AR) w 
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For simplicity, confine ourselves by the case / = . 
From relation ^ whose both sides in this case con- 
cern equilibrium (hence, spherically symmetric) Brow- 
nian motion, we can connect V\ (t, AR) — V\ (t, AR ; 1) 
with derivative of the left side by the relative gas density: 
Vx(t,AR) = [ dVo (t, AR ; 1 + 0) fd<j> ] ^ = o . Then, com- 
bining this equality with @ and noticing that functions 
F\ and Wi by their birth from DF F\ are non-negative, 
in view of Wi > we come to inequality 



Cl V {t,AR;l) 



dV (t,AR;l + <P) 



> (10) 



> = o 



Further, let us assume that at t 3> t = A/vq , 
where r is mean free flight time and Vq characteris- 
tic thermal velocity of BP, a probability distribution of 
BP's path tends to the Gaussian one (see paragraph 1): 
Vo(t, AR) — > Va(t,AR) . Besides, take into acount that 
diffusivity of BP, D , in the Boltzmann-Grad gas is ar- 
ranged as D — Avq = VqT and, consequently, together 
with A is inversely proportional to gas density. There- 
fore left side of ((U should follow from Va{t,AR) after 
change A -» A/(l + <f>) , D — > D/(l + (f>) , that is 



V (t,AR;l 



1 



4irDt 



\ 3/2 

J exp 



-(1 



AR 2 

ADt 



Insertion of this expression to inequality (fTU|) yields 
V G (t,AR) [d + 3/2 - AR 2 /4Dt] > , i.e. produces 
obvious contradiction. 

Hence, the theory does not allow for the Gaussian 
asymptotic! This statement, or rather its new proof since 
it itself is not novell (see paragraph 1), is principal result 
of the present paper. 

Of course, it is not rigorous because we have managed 
without a concrete definition of the averaging procedure 
marked by the over-line (though hcuristically this is ad- 
vantage of our consideration). But, in view of ([7]), this 
procedure may influence upon numerical value of the co- 
efficient ci . This notion, nevertheless, does not spoil our 
result, since it means merely that we should substitute to 
pop minimum of possible values. Our presumption that 
ci is indeed a coefficient, but not a function of t, AR, Po 
and pi in ([7]) or t and AR in ([5]), also can be justi- 
fied. Physically, the length of collision cylinder, ciA , 
for pair of particles is bounded by collisions with "third 
particles" which knock out the pair from a number of 
candidates for collisions. Since the "third particles" orig- 
inate mainly from uniform background, the constancy of 
Ci seems quite natural. 



5. In conclusion, a few more remarks. 

First, the above formulation of the Boltzmann-Grad 
limit implied that size of BP is comparable with size of 
atoms. Therefore conclusions of previous paragraph do 
not apply to "macroscopic" BP (in respect to it, at least 
in the limit of infinitely large BP's mass, Gausian asymp- 
totic remains above suspicion). 

Second, and what is allowed for by the inequality 
(fTOf ? Let Vo(t, AR) at t > r is characterized, simi- 
larly to Va(t, AR) , by a single parameter, that is dif- 
fusivity: V (t,AR) -> (2Dt)' 3 / 2 ^(AR 2 /2Dt) , where 
/ ty(a 2 ) da = 1 . Then inequality (fT0|) reduces to 

(ci +3/2)*(x) + xd^{x)/dx > 0, 

thus requiring that function ^(x) must have power- 
law long tail: ^(x) cx x~ a at x — > oo , 
where, clearly, a < c\ + 3/2 . But such behav- 
ior would lead to unboundedness of higher-order sta- 
tistical moments of BP's path. Recall, however, that 
we have once more parameter, vq , and hence, at 
t » r, the possibility to write Vo(t,AR) — > 
(2Dt)- 3 / 2 $>{AR 2 /2Dt)Q(\AR\/v t) , where 6(0) = 1 
and O(x) enough rapidly decreases at infinity. From the 
point of view of (JTUJ) , this is also allowed variant. And it 
reproduces, if take c\ = 2 , the asymptotic of distribution 
of BP's path found in [12| in the frame of "collisional ap- 
proximation" to BBGKY equations [Hi EH- At that, the 
role of BP was played by test (marked) atom. Such the 
asymptotic reflects non-ergodicity of its Brownian tra- 
jectories, which can be characterized as "flicker" fluctu- 
ations of diffusivity (and mobility) of BP [ljj [H H3 . 

It should be underlined that our present formalism 
and one exploited in [l2j have nothing common among 
them (except their application). Therefore the remark- 
able nearness of their results can be rated as evidence for 
adequacy of both approaches. 

At last, thirdly, in our above consideration by the ex- 
ample of pair DF we in fact have outlined connections 
between the correlation functions, on one hand, and the 
specific distribution functions (for pair and many-particle 
collisional configurations) first introduced in [1JJ] , on the 
other hand. In particular, function F\ here is equivalent 
of function A 2 from [l(| , while function W\ defined in 
© is equivalent of W% from pjj]. Extension of these 
connections to higher-order DF will unify the two for- 
malisms and, undoubtedly, open new applications of the 
virial expansions ((4| and (J6j> . 

I would like to acknowledge Dr. I. Krasnyuk for many 
useful conversations. 
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